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O Abstract 

A general class of solutions is obtained which describe a spherically 
symmetric wormhole system. The presence of arbitrary functions al- 

q -1 lows one to describe infinitely many wormhole systems of this type. 

^ \ The source of the stress-energy supporting the structure consists of 

PP' an anisotropic brown dwarf "star" which smoothly joins the vacuum 

and may possess an arbitrary cosmological constant. It is demon- 
k>( \ strated how this set of solutions allows for a non-zero energy density 

I— , ' and therefore allows positive stellar mass as well as how violations 

of energy conditions may be minimized. Unlike examples considered 
thus far, emphasis here is placed on construction by manipulating the 
matter field as opposed to the metric. This scheme is generally more 
physical than the purely geometric method. Finally, explicit exam- 
ples are constructed including an example which demonstrates how 
multiple closed universes may be connected by such wormholes. The 
number of connected universes may be finite or infinite. 
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1 Introduction 

The study of spacetimes with nontrivial topology is one with an interest- 
ing history. One of the earliest, and most famous, wormhole type geometries 
considered is that of the Einstein- Rosen bridge |l| . This "bridge" connecting 
two sheets was to represent an elementary particle. That is, particles both 
charged and uncharged were to be modeled as wormholes. From the work of 
Einstein and Rosen one also can find a hint to what eventually was to become 
a central problem in wormhole physics. Energy conditions must generally be 
violated to support static wormhole geometries. 

The wormhole resurfaced later in the study of quantum gravity in 
the context of the "spacetime foam" |2|] || in which wormhole type struc- 
tures permeate throughout spacetime at scales near the Planck length. More 
recently, there is the meticulous study performed by Morris and Thorne |4| 
considering a static, spherically symmetric configuration connecting two uni- 
verses. Since the work by Morris and Thorne, there has been a considerable 
amount of work on the subject of wormhole physics. For an excellent review 
the reader is referred to the thorough book by Visser as well as the exhaus- 
tive list of references therein. Most interesting is the case of the traversable 
wormhole which, at the very least, requires tidal forces to be small along 
time-like world lines as well as the absence of event horizons. These geome- 
tries are of interest not only for pedagogical reasons but they serve to shed 
light on some fundamental issues involving chronology protection || and the 
topology of the universe. Also, many wormhole solutions are in violation of 
the weak energy condition and there has been much activity in the literature 
to attempt to minimize such violations (see, for example, [0, [§], §|, flO 



and references therein.) There is also the possibility that wormholes have 
much relevance to studies of black holes. For example, it may be possible 
that small wormholes may induce large shifts in a black hole's event horizon 
or that semiclassical effects near a singularity may cause a wormhole to form 
avoiding a singularity altogether jn{] , JT2] , ]T5| . 

For the above reasons, and many others which are hinted at in the 
literature, it is instructive to study as general a class as possible of such solu- 
tions. The aim here is to construct a mathematical model which encompasses 
the bulk of static, spherically symmetric, wormhole solutions. This is done by 
studying the profile curve of a generic wormhole structure and postulating a 
mathematical expression describing this curve. From this ansatz all relevant 
properties may be calculated. Here we are concerned with making the system 



as physically relevant, yet as general, as possible. This involves, as much as 
possible, working with the stress energy tensor when studying the system, 
as opposed to fully prescribing the geometry as is usually done. We demand 
that the physical requirements be as reasonable as possible given the situa- 
tion and study the system from this point of view. Also, the system is not 
supported by a matter field which permeates all space but possesses a definite 
matter boundary. At this boundary the solution is smoothly patched to the 
vacuum with arbitrary cosmological constant. The presence of a cosmologi- 
cal constant is permitted for several reasons. First, astronomical studies of 
distant supernovae favour the presence of a non-zero cosmological constant 



I 1] p!5| . Second, the cosmological constant necessarily violates the strong 
energy condition and therefore may serve to minimize energy condition vi- 
olations of the matter field. We make comments where appropriate on how 
the introduction of the cosmological term affects the behaviour of the matter 
field. In this paper we find that solutions may be obtained in which the 
matter field does not violate energy conditions at the throat and with min- 
imal energy condition violation elsewhere. It is possible that a patch to an 
intermediate layer of matter may eliminate violations altogether. 

The paper is organized as follows: Section 2 introduces the Einstein 
field equations for a general static spherically symmetric system, as well 
as the differential identities which must be satisfied. This system is then 
specialized to describe a general wormhole system. Section 3 illustrates the 
methods used to solve the field equations. The matter field is also introduced 
and restrictions on its behaviour are analyzed. The singularity and event 
horizon conditions are studied here in terms of the matter field. The smooth 
patching of this solution to the exterior Schwarzschild-(anti) de Sitter solution 
fl~6| is also covered here along with junction conditions. Specific examples 
are constructed and energy conditions are studied in section 4. Finally, some 
concluding remarks are made in section 5. 

2 Geometry and Topology 

The geometry studied is displayed in figure 1. For the moment we restrict 
our analysis to the case where there exist two regions, an upper and lower 
sheet, which are connected by a throat of radius r = tq. Neither region 
need be asymptotically flat and each region may represent a universe which 
is connected by the wormhole throat. It would not be difficult to modify 



the calculations here to represent a throat connecting two parts of the same 
universe. This, for example, may easily be accomplished via topological 
identification at spatial infinity in the case of an open universe. The closed 
universe case will be considered as a specific example later. 

In the case studied here both the upper and lower universes possess 
a matter-vacuum boundary at r = d± > r where the + applies to the 
"upper universe" and the — denotes quantities in the "lower universe" . This 
boundary represents the junction between the brown dwarf and the vacuum. 
We use the term brown dwarf to describe the matter system for several 
reasons. First, the traversability condition precludes the use of a matter 
system undergoing intense nuclear processes in its interior. A very rough 
definition of a brown dwarf is a "star" which forms via gravitational collapse 
but does not ignite nuclear fusion in its core. Second, it is reasonable that if 
such a system were to form there would be some cut-off point in the stress- 
energy Our system here possesses such a boundary (the stellar surface) and 
therefore the analogy with stellar structure again leads us to the term brown 
dwarf. It should be noted that d + does not necessarily have to be equal to 
gL and therefore the brown dwarf "star" may have a different size in each 
universe. 

We take as our starting point the spherically symmetric line element 
of the form 

d s 2 = -e 7(r)± dt 2 + e a{r)± dr 2 + r 2 d6 2 + r 2 sm 2 (9) dtf (1) 

with j(r)± and ot(r)± being functions of the radial coordinate only. The 
coordinates cover the range 

-oo < t < oo, < r < r < oo, < 9 < n, < < 2n. (2) 

This form will prove especially useful when patching solutions to the vacuum 
at r — d±. The condition that there be no horizons implies that j(r)± be 



bounded |T7|]. Also, we demand that the time-like coordinate, t, be smooth 
across the throat which yields the condition 7(r )_ = 7(r ) + . 

The fundamental equations governing the geometry are the Einstein 
field equations [] with cosmological constant 

RtMu - -Rg^v + Afl> = 871-T^, (3) 



1 Conventions here follow those of j| with G = c = 1. The Riemann tensor is therefore 
; iven by W pva = r^„ + T^T a pa -... with R pa = R a paa . 
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Figure 1: Embedding diagram for the wormhole geometry separating two uni- 
verses. Circles of constant r represent two-spheres. There is a throat at r = tq 
and the boundary of the brown dwarf is located at r = d± where it is patched 
smoothly to Schwarzschild-(anti) de Sitter geometry. 



along with the conservation law 

T^ u = 0. (4) 

The spherically symmetric Einstein field equations, in mixed form, 



yield only three independent equations: 

Gl = 6 -^- (1 - ra(r) ±il ) - 1 = SnT ± - A (5a) 

Gi 1 =^^ (1 + r T (r) ±jl ) - 1 = SttT^ x - A (5b) 

G 2 ± 2 = G| 3 =^^ ( 7 (r) ± ,i,i + ~ (7(r) ± ,i) 2 + J (l(r)± - <*{r) ± ) tl 

ia(r)±,i 7 (0±,i J = 8ttT1 2 - A = SttTJ! 3 - A, (5c) 



2 
whereas the conservation law yields one non-trivial equation: 

^1,1 + (l + j7(r),i) ri - ~7(r),iT° = ^ 2 = ^3- (6) 

Note that the cosmological constant, A, is assumed to have the same value 
in both the + and — region. This is because the cosmological constant is 
taken to represent the stress-energy of the vacuum and we assume here that 
both regions, being connected, possess similar vacua. It may be interesting 
to study situations where, for example, one universe is of de Sitter type and 
the other of anti-de Sitter type. In this case the de Sitter universe would be 
a "baby universe" to the anti-de Sitter one. 

Equation (Sal) may be utilized to give the following: 
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- ^Vffc -A),^ + 1 . :1 .^ (7) 

r J ro V 87r / r r 

with b(r)± defined as the shape function. If one is considering standard 
stellar models with plain topology, one must set the constant c = to avoid 
a singularity at the origin. Here, however, the origin (r = 0) is not a part of 
the manifold (see figure 1) and this constant must be set by the requirement 
that the throat region be sufficiently smooth, as described next. 

Figure 2 shows a cross-section or profile curve of a two dimensional 
section, t = constant, 6 = n/2, of the wormhole near the throat. The 
wormhole itself is constructed via the creation of a surface of revolution when 



one rotates this figure about the x 3 -axis. The surface of revolution may be 
parameterized as follows 

x(r,(f>) := (x x (r, 0), x 2 (r, 0), x 3 (r, 0)) = (rcos(0), rsin(0), P(r)) , (8) 



where r = \/ (x 1 ) + (x 2 ) and is shown in figure 2. Therefore, the induced 



metric on the surface is given by 

da 2 = [l + Px(r) 2 ] dr 2 + 



,2 ji2 



(9) 



where commas denote ordinary derivatives. Note that the corresponding 
four-dimensional spacetime metric, ([!]), may now be written as: 



ds 2 = -e 7(r) dt 2 + [1 + Pi(r) 2 ] dr 2 + r 2 dd 2 + r 2 sin 2 (^) 



(10) 



The shape function, b(r), may now be specified in terms of the embedding 
function, P(r), which will be discussed shortly. 



x 3 = P(r) 




Figure 2: A cross section of the wormhole profile curve near the throat region. 
The wormhole is generated by rotating the curve about the x 3 -axis. 



To set the constant, c, in (0) note that from the profile curve in 
figure 2 the derivative Pi(r) tends to infinity as one approaches r . This 
condition translates, via ([TOD, to 1 — b(r)±/r vanishing at r = tq. Therefore 
the constant c must be equal to — vq. This condition is equivalent to that 
in [|J and reflects the fact that at the throat a finite change in coordinate 
distance, Ar, corresponds to an infinite change in proper radial distance, Al. 
That is, 

Ar 
lim(AZ) = lim — ^=^= -> oo. (11) 

*ro _ / 1 _ b{r) ± 



A final condition for smooth patching is the continuity of the second funda- 



mental form ||18| . For the moment we assume that the joining of the upper 
and lower universes may be performed. It will be shown, in section 3, that 
the second fundamental form is indeed continuous however; a few more prop- 
erties of the system need to be studied in order to demonstrate this. 

We are now in a position to analyze the geometry. For now we limit 
the analysis to the top half of the profile curve (the "+" region) since the 
bottom half is easily obtained once we have the top half. From the profile 
curve in figure 2 it may be seen that a function is needed with the following 
properties: 

• The derivative, -Pi(r), must be infinite at the throat and finite else- 
where. 

• The derivative, -Pi(r), must be positive at least near the throat. 

• The function P(r) must possess a negative second derivative at least 
near the throat region. 

• Since there will be a cut-off at r — d, where the solution is joined to 
vacuum (Schwarzschild-(anti) de Sitter), no specification needs to be 
made as r tends to infinity. 

These restrictions are purely geometric in nature. We show below that rea- 
sonable physics will demand further restrictions. From these properties a 



general function is assumed of the form: 

P(r) = A(r-r ) e exp[-/i(r)], (12a) 

P tl (r) = P(r) [e (r - r^ 1 - h,i(rj\ , (12b) 

P(r) xl = P(r) { [e (r - r^ 1 - h tl (r)] * - h(r) xl - e (r - r )- 2 } . (12c) 

Here h(r) is any non-negative (at least twice differentiable) function with 
non-negative first derivative. Moreover, e is a constant bounded as follows: 
< e < |. The reason for this restriction will be made clear later. A is a 
positive constant which is required so that quantities possess correct units 
(the unit of A is length 1-6 .) One may enforce conditions on h(r) by writing 
it in terms of a slack function as 



h,i(r) = e 



q(r) 



r 



h(r) = / e q{r>) dr' + z , 



'■() 



with zq a constant and q(r) an arbitrary differentiable function. 
The shape function, b(r), is given by (0) and fllOD 



& (0 = ^/p V2 , (13a) 

1 + P A (ry 
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P,(r)+ 2rP(r) ^ 



1 + Pi(r) 2 



(13b) 



The bottom universe may be obtained via the profile curve by demanding 
that A be some negative constant. The first derivative will then be negative 
near the throat and the second derivative positive. The function h(r) may 
be different in the lower universe than in the upper and this will not spoil 
smoothness at the throat as long as h + (r ) = /i-(r ) and /i +i i(r ) = 7t_ i(ro). 
To show that the system is indeed a wormhole with a smooth patch- 
ing at r , it is useful to show that the system may be described by a single 
coordinate chart across the throat. Consider the profile curve in the x l x 3 - 
plane of figure 2. The upper half of this curve may be described by 

x 1 =X\ (r+) := r+, (14a) 

x 3 =X% (r + ) := P+(r+), (14b) 



with the + subscripts indicating the upper universe with r + > r . We denote 
the inverse of P + {r + ) from ( |14b| ) as r + = Pl(x 3 ) (Jj superscripts will be used 



to indicate inverse mappings). Also, the lower half of the curve is given by 

x 1 =Xi(r_) :=r_, (15a) 

x 3 =Xi(r_) :=P_(r_), (15b) 

and again r_ > r® with r_ = P£(x 3 ). At the moment we have two charts 
in two universes and we may construct two metrics from these (the "+" and 
"-" metrics from above). The unifying chart, which includes both universes 
near the throat as well as the boundary, tq, is given by studying the profile 
curve as a function of x 3 . In other words, we rotate the x 1 x 3 -plane about the 
x 2 -axis of figure 2 by tt/2 and study the profile curve parameterized by the 
coordinate x 3 . By equation (|12a|) , x 3 belongs to an interval containing x 3 = 
which corresponds to the throat of the wormhole (r corresponds to P^(0)). 
With this parameterization, P"(0) = r$ exists as well as lim a .3^ 0+ P\x 3 ) = 
lim^^o- P^(x 3 ) and therefore P" is continuous. The functions P|(x 3 ) and 
P_(x 3 ) are simply obtained by the restrictions of P"(x 3 ) in x 3 > and x 3 < 
respectively. The resulting curve is parameterized as: 



x 



X*\x*) = P\x% (16a) 



x 



X i3 (x 3 ):=x 3 . (16b) 



With this "universal" chart the spacetime metric near the throat may be 
written as 



ds 2 = _ e Y^) dt 2 + [ 1+ py^ ( rfx 3)2 + p8 (x 3 ) 2^2 + ptt (x 3 ) 2 sin2 ^ 2_ 

(17) 



,x" 



This single metric is therefore sufficient to describe both parts of the worm- 
hole system across the throat. 

In the coordinate chart used thorughout most of this paper, (|TUD, 
the manifold must be at least piece-wise C 3 smooth to have a well behaved 
Einstein tensor with contracted Bianchi identity. At the throat junction 
we will later show that a minimum of C l is admitted although there is no 
reason why a higher, perhaps more physical, class may not be admitted. 
It may be argued that the chart given by ( pf ) is better suited for studies 
across the throat. In that case the embedding function, P"(x 3 ), may be 

10 



chosen of arbitrary differentiability class (ideally minimum C 4 to give a C 3 
manifold) and no piece- wise considerations need be addressed at the throat. 
The coordinate transformations from one chart to the other, however, may 
turn out to be formidable. The above discussion refers to the differentiability 
of the metric. Since the metric depends on the derivative of the profile curve, 
a C n manifold should possess a C n+1 profile curve. 

3 Wormhole Structure 

In this section the field equations are utilized to construct the wormhole 
structure. There are a number of ways to proceed and each method has its 
particular advantages and disadvantages. One method is a purely geometric 
method or g- method [[HJ. In this method one is free to specify the metric 
components, j(r) and b(r), to acquire the desired geometry. One then com- 
pletely determines the supporting stress-energy tensor from the geometry via 
the field equations. This is the method of Q. One may, for example, wish 
that j(r) be bounded to eliminate horizons. As a simple example, such a 
function may be expressed in terms of a sum of an even and an odd bounded 
function as: 

7(7") = c\ tanh(</?(r)) + c^ sech(^(r)), (18) 

which is bounded as 

| 7 (r)| < |ci| + |ca|. (19) 

Here tp(r) and ip(r) are arbitrary. This method possesses the advantage that 
any reasonable geometry allowed by general relativity may be constructed. 
The disadvantage, however, is that one may be left with a stress-energy tensor 
which is physically unreasonable. 



A second method is the T-method |I9 |. Here one prescribes proper- 
ties of a desired matter field and from this constructs a stress-energy tensor. 
The field equations are then solved for the metric which governs the geom- 
etry. This is generally a preferred method as one may prescribe physically 
reasonable matter and study the resulting gravitational field. The disadvan- 
tage, of course, is that there is very little control over the geometry. If we 
have a particular geometry in mind, there is no guarantee that an a priori 
prescribed stress-energy tensor will generate the desired gravitational field. 
Therefore, for the purposes of wormhole analysis, this method is not useful. 

11 



A third method, and one which most of the analysis here is based 



on, is a mixed method [20]. This method may be implemented in two ways: 



• Prescribe some desirable physical properties to T_j? and T±±. For in- 
stance, for physically reasonable matter one may wish the first to be 
negative and the latter non-negative. Otherwise, T± ± may be freely set. 
T± may be obtained via the shape function, b(r)±, using (|5a]). The 
physical properties demanded above must be utilized to constrain the 
shape function. The transverse pressures T± 2 = T± 3 are defined via (||). 

• Again obtain T± from the shape function as above and prescribe an 
equation of state relating T_y to T± Q . Again the transverse pressures 
are obtained from (|6]). 

All equations will then be satisfied. The advantage of this method is that 
one may employ some reasonable physical assumptions on the matter fields, 
yet not give up all control over the geometry. The major disadvantage is 
that, since g 00 is not set, certain desirable properties of the manifold, such 
as the absence of horizons and curvature singularities, are not guaranteed. 
These problems will be addressed below. 

The matter field may, for example, be that of an anisotropic fluid 
whose stress-energy tensor is given by 

Tp, = (p + p±) u^u u + p ± g^ + (p„ - p±) s^s v (20) 

with 

u» Ull = - 1, (21a) 

s% =1, (21b) 

w% =0. (21c) 

Here, in the static case, the eigenvalues of T£ are T° Q = —p, T\ = p u and 
T\ = T\ = p± which are a measure of the (negative) energy density, radial 
pressure and transverse pressure, respectively. This matter model fits our 
prescription for several reasons. First, the construction of static, spherically 
symmetric wormholes requires a stress-energy tensor which admits four real 
eigenvalues. Two of these eigenvalues coincide by demanding spherical sym- 
metry [T\ and T|). Second, utilizing both the g and mixed methods will most 
likely result in a system possessing three distinct eigenvalues. Furthermore, 
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it is known that these eigenvalues must, at least in some vicinity, violate 



energy conditions |4j]. It has been shown |21[ that the anisotropic fluid, when 
undergoing gravitational collapse, may form many states of "exotic" (energy 
condition violating) matter. Finally, a wormhole system as described in this 
paper would most likely form via gravitational collapse, as with normal stars. 
The anisotropic fluid therefore represents a matter field which is a reasonable 



extension of certain standard stellar models involving a perfect fluid f22 |. 
The metric function , j±(r), is given by (|5b|) and its integral: 



TW±,i 
7(0± 



SttT^W-A 



r 



£H iM - A+ w)i 



i 

b±(r) ~ r' 

r' dr' 



b±(r>) 



(22a) 
ln(r) + K±, (22b) 



with K± a constant. Continuity of j(r)± at the throat requires that K^ 
K_; a(r)± is given by equation (|7|): 



a(r)± 



-In 



1-^ + ^ 



J\2 



''() 



rpO 

1 ±0 



A 



57T 



dr' 



a(r) 



In I 1 

1 



6±(r) 



±.i 



Mr) 



K r )±,l - 



b{r) 



(23a) 
(23b) 



Another useful equation, which will be utilized later, is obtained by subtract- 
ing fl5al) from (|5b|) 



Svre^r (T^ - T° ) = ( 7 (r) ± + a(r) ± ), . 



(24) 



Having expressions for the metric functions and their derivatives, we 
now turn to the matter field. This analysis will allow us to determine all 
properties which the general embedding function need to satisfy and will 
therefore provide a general class of metrics describing such wormholes. 

The Einstein equation fl5a|) may be written in terms of the shape 
function yielding: 



1 



7 b{r), 1 = torT£ -A. 



(25) 
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We assume here that, at least in the vicinity of the throat where the cos- 
mological constant is negligible compared to the matter energy density, the 
net energy density is positive. From ( f25| ) it is seen that this requires that 
b(r) } i > near the throat. Recall from the properties of the profile curve 
(figure 2), P(r) must possess positive first derivative and negative second 
derivative near the throat. From ( |13b| ) it may be seen that if the second 
term in the expression dominates, b(r) t i will be negative which is not desir- 
able. Therefore (|13b| ) must be analyzed with some care in the vicinity of the 
throat. 

To study the properties of ( p5[ ) near the throat, we define a parameter 
5 := (r — r ). Analysis of &(V) i then shows that 

b{r) tl nl-0{6 1 - 2e ) + ... (26) 

from which it may be seen that e must be restricted to the interval < e < | , 
as mentioned earlier. Note that e = ~ is allowed, although a negative energy 
density as observed by static observers is the price one pays for this case. 
Even for the values of e stated above &(r) i eventually becomes negative at 
some radius, r = r± say, as is evident from ( |13b| ) . One could, if one insists on 
a positive value of b(r) t i everywhere, place the stellar boundary at d± < r±. 
Note that since the model admits an arbitrary cosmo logical constant, we 
could allow b(d±) i i < and still claim a positive p for the matter field if 
A<0. 



The equation (p6|) with the above restriction on e yields the best 
possible scenario as far as the energy density is concerned. That is, a value 
of b(r)± t i greater that unity cannot be achieved at the throat. This situation 
is quite favourable with respect to the energy conditions (to be discussed 
section 4) as the null, weak and dominant energy conditions each require 
positivity of energy density. 

3.1 Horizons and Singularities 

As mentioned earlier, the method utilized here, although giving some 
control over the matter fields, does not guarantee the absence of event hori- 
zons nor singularities. In this section we address both of these issues. The 
notation r(r) := 8ttT\ — A is used throughout this section. 

In the coordinate system of ([I]) an event horizon will exist wherever 
goo vanishes. In other words, the function 7±(r) must be everywhere finite. 
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Using ( j22a|) and ( j22b| ) it may be seen that, to ensure that r y{r)± remain finite 
near the throat, (Y(r) + 1/r 2 ) must vanish at least as fast as (r — b(r)±) in 
the limit r — > tq. 



From (|13a| ) it may be shown that near the throat 

(6(r) - r) oc (r - r ) 2(1 " e) ; (27) 



therefore the numerator of the first term in (|22a|) must possess behaviour 



near the throat of the form (r — r^)" 1 with r\ > 2(1 — e) and 

r(r ) = -l/r 2 , (28) 

Although we may possess everywhere non-negative energy density, (^) shows 
that a negative pressure or tension is required at the throat to support the 
wormhole. The fluid pressure may still be positive if the cosmological con- 
stant is positive and dominates the fluid pressure at the throat. 

The equation ( |22b|) may most easily be made well behaved for any e 
allowed by the model by demanding that the first-order term in the expansion 
of r(r) + 1/r 2 about r$ vanish. Also, the appropriate limit r(ro) = — I/Vq 
must be enforced so that the zeroth-order term vanish as well. With little 
loss of generality we may write r(r) as 



r(r) = — - cos 



tq) rmr 



cosh [(r — r ) s £(r)] 



r 2 l(d + a — r ) 2 
+ K(r-r ) q (r-d) w exp(((r)) (29) 

with m an odd integer, q > 2,w > positive integers and £(r) and ((r) 
arbitrary functions save for the restriction that they possesses well behaved 
derivatives in the regime r$ < r < d; a is a small constant set by the 
requirement that at r = d± the above expression is equal to —A. The second 
term, proportional to the constant k, is not strictly required, though it is 
useful in making r(r) positive throughout most of the stellar region (see 
section 4). An analysis of the matter field in the context of the energy 
conditions may also be found in section 4. 

Having eliminated the possibility of event horizons, it is now neces- 
sary to analyze the manifold for singularities. Singular manifold structure is 
most easily studied by constructing the Riemann tensor in an orthonormal 
frame. We use notation such that indices surrounded by parentheses denote 
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quantities in the orthonormal frame. In this frame the Riemann tensor pos- 
sesses the following components as well as those related by symmetry 0: 



R, 



(trtr)± 



6±(r) 



l{r)±,i,i + -(j{r)± >1 ) 2 



1 ( , 6(r) ± ,i - ^ 



i2(tetfl)± =2 _ 7('")±,i I ! - 



r — 6±(r) 
6±(r 



R 



1 
6±(r 



6(r)±,i - 



r 

6±(r) 
r 



with 



- R(t<j)t<t>)± 



and 



6(r ) 



(30a) 

(30b) 

(30c) 

(30d) 
r 



Recall that, utilizing our ansatz 
and b(ro),i = 1. Therefore all components of the Riemann tensor, with 
the possible exception of ( |30a| ), are finite at the throat. The only other 
problematic spot is r = 0, which is not a part of the manifold and therefore 
causes no concern for our analysis. The previous restrictions on r(r)± ensure 
that j(r)± and its derivatives are finite away from the throat. No quantity 
grows without bound as we move away from the throat. 

To study ( p0a| ) more carefully, it is useful to write it in terms of T^ 
using (|22a|) as (dropping the ± subscript for now) 



R 



(trtr) 



4 r 3 ( r _ b(r)) 



[r(r) 2 r 6 + 2r 4 r(r),i (r - b{r)) 



_l_ r 3 r ( r ) (2r — b(r) + rb{r)i) — Arb{r) + 2r 2 5(r) 1 
+ 46(r) 2 - b(r)b(r),ir] . 



(31) 



Notice that for non-singular behaviour, the term in square brackets must 
vanish at least as fast as the denominator as r — > r$. Recall that the denom- 
inator's behaviour is as in ( p7|) and therefore the numerator must possess 
similar or stronger vanishing properties here. An exhaustive calculation us- 
ing (|29"D reveals that the numerator does indeed vanish at least as fast as 
the denominator near the throat and therefore all curvature tensor compo- 
nents are finite, eliminating any singularities. It is interesting to note that, 
by demanding that j(r)± t i be finite everywhere (the horizon condition), we 
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also impose the final condition which ensures finiteness on all orthonormal 
Riemann tensor components. In other words, naked singularities are not 
supported. 

We now return to the consideration regarding the continuity of the 
second fundamental form at the throat. The second form is most easily 
calculated in the coordinate frame as the covariant derivative 

Kin, = - {n KV + hu.fjt) , (32) 



where n M represents an outward pointing radial unit vector which, using (|ID| 
is given by: 



n, 



1 + P,i(r) 2 *W (33) 



p ~ V ' > lv * a* 



(the ± subscripts have been dropped here). From (|3^) the following compo- 
nents are calculated for the second fundamental form: 

K tt = - —= L= e^ 7 ,i(r) , (34a) 



2^/1 + P tl (r 



|2 



K ^^WW' (34b) 



rsin 2 ^ 

K » = TTWTT 2- (34c) 

1 + PAry 



( 34b ) and (|34cj ) are certainly continuous at the throat since, from ( |12b|) they 



both vanish as one approaches r = r from either the "+" region or the "-" 
region. K tt will be studied with some care. Recall that the condition ensuring 
the absence of horizons is given by demanding that (|22a|) be everywhere finite. 



This condition therefore requires that t{tq) := 8irT\(r) —A = — 1/ro and this 
result is, of course, independent of the region in which one is approaching tq. 
This condition, along with the fact that b + (ro) = 6„(r ) = r dictates, via 
( 22a ), that 7,1(7") be continuous at the throat. Also, continuity of the metric 



requires that e 7 ^ r °' )+ = e 7 ^ r °^ and therefore K u must also be continuous at 
the throat. All required junction conditions at the throat have therefore been 
met (namely: continuity of the metric, T\(r), and the second fundamental 
form as well as condition (|ll]) along with visual proof from figure 4b to be 
discussed in section 4). 
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If one is considering a traversable type wormhole then a further re- 
striction on Riemann components is generally desired. Namely, the magni- 
tude of these components must be small in the sense that tidal acceleration, 
given by: 

^ = K 0J V a V p AC (35) 

be sufficiently small to allow travel through the throat. Here A£ 7 is a sep- 
aration vector joining two parts of the traveling body. The condition that 



( 55]) be small places constraints on the magnitude of the orthonormal Rie- 
mann tensor components as well as on the velocity at which the wormhole 
may safely be traversed 0, |5[]. At the throat, it is possible to make all 
components of (|30a|) -( J30dD small or vanish due to the fact that b(r ) = r 
and b(r )i = 1, with possible exception of ( ^0d| ) which is equal to I/tq. For 
purely radial motion ( |35"D dictates that this component plays no part in the 
tidal acceleration and therefore presents no constraint on the wormhole size 
H, Pj. Away from the throat there is the potential to make the relevant 



components at least reasonably small by the choice of matter field via (|22a|) 
and fl3T|). Note that the "smallness" of the Riemann components is not a 
requirement of the wormhole system. It is simply a measure of the ease at 
which one may be traversed. 

3.2 The Stellar Boundary 

Having identified conditions in the interior we now turn attention to the 
boundary of the star at r = d±. Here the solution must smoothly tend to 
the vacuum Schwarzschild-(anti) de Sitter metric: 

2M 0± 1 A ,_ dr 2 

ds z = - 1 — --Ar 2 )dP + 



3 J ^_2A^_I Ar 2 

+ r 2 d6 2 + r 2 sin 2 (0) d<f) 2 (36) 

where the tilde on the time coordinate will be later made clear. Note that 
the "mass" of the star in one universe need not necessarily be the same as 
in the other. The junction condition employed at this surface will be that of 



Synge [19 



nX\ r=d+ = (37) 



with n v an outward pointing radial normal vector. This condition essentially 
states that there may be no flux of stress-energy off the surface of the star 
which, of course, is a reasonable requirement when patching to a vacuum 
solution as will be done next. For the class of metrics considered here, con- 
tinuity of the metric at the boundary also ensures continuity of the extrinsic 
curvature. 

It is now necessary to check conditions at the boundary r — d± where 
the interior solution must match the vacuum solution given by d36|). Here, 
from (0) we get 

-l 



P a(<2±) 



1 



cL 



A r% 
3 d± 



2M± + a r$_ _!±_tl d 2 



(L 



A 
— i 
3 



where M± corresponds to the fluid mass 

M± = -4vr 



r'fTUdr'. 



(38) 



(39) 



The function ( pB]) matches smoothly to the Schwarzschild-(anti) de Sitter 
metric (B6|) with an "effective mass" 



Mr 



o± 



M± - '-r 3 



A 
6" 



To 

2 ' 



(40) 



The equation for r j± is more complicated: 

2M 0± A 



,7(4t) 



x exp 




(r') 2 (p±„ + p±) dr' 



'-87r/:(r*) 2 (p± + ^)dr*-r 



1 



2M 



o± 



A 



J: 



3 2x± 



+ K 



4= s-J- - (41) 

Here (p4]) has been used. We now perform the following coordinate change 

i = e x t (42) 

so that 

(di) 2 = e 2x dt 2 (43) 

and thus the interior solution has been matched to (1361). 
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4 Examples 



Having developed the formalism for constructing such wormholes it is in- 
structive to demonstrate how the scheme works via specific examples. Again 
the ± subscript is neglected unless needed for clarity. 

A simple, yet instructive example is one with the following parame- 
ters: 

h{r) oc r 2 , 

£(r) oc r 3 , 

m = 3, 

8= 1, 

A ~ small and positive, 
K ~ large and positive, 
((r) a constant. 

The rotated profile curve for this system is shown in figure 3a. Figure 3b 
displays the shape function and 3c its derivative. Note that the total mass 
(as measured by b(r)) and the energy density (&(?"), i) are positive in the 
vicinity of the throat. Since the Einstein equations permit discontinuities in 
the energy density, there is no reason to place the stellar boundary near the 
point where b(r) t i vanishes although for this example we do so. The radial 
pressure, p n = (r(r) + A) /8n, is plotted in figure 3d. Note that with the 
addition of the large k term this term is positive throughout most of the 
region. Finally, figure 3e shows a plot of the transverse pressure. As with 
the radial pressure, this quantity is positive throughout most of the stellar 
bulk. 



4.1 Energy Conditions 

Here we discuss energy conditions both in general and with respect to 
the above example. From an expansion of (E5l) near the throat we find that 






A 2 e 2 r 2 
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(3c): db(r)/dr 



(3b): b(r) 



(3d): radial pressure 



(3e): transverse pressure (3f): weak energy coll dition 



Figure 3: A sample wormhole system, a) shows a rotated profile curve of the "+" 
universe, b) is a plot of the effective mass of the star, c) a plot of the effective 
energy density, b(r) t \, d) the parallel pressure e) the transverse pressure and f) the 
weak energy condition. Note that p is non-negative. Moreover, throughout most 
of the stellar bulk p u and p± are positive. 



A similar analysis on the radial pressure yields via (p9|): 



Sirp n 



+ A + ^(r-r ) + O(r-r ) 



(46) 



Finally, the transverse pressure 
87rp_Lv.ro) 



^* o Pi M \r=r 



A = A. 



(47) 



Exactly at the throat, energy conditions are respected as p is positive unless 
the cosmological constant is large (a possibility ruled out by experiment) 
and p(r Q ) + p u (r ) = as well as p(ro) +p±(r ) > 0. As we move away from 
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the throat notice that, in the expression p + p u , the third and fifth terms of 
(|45|) survive. Both of these terms yield a negative contribution and therefore 
energy conditions are violated slightly away from the throat. It is possible 
to cut-off the solution at some radius near r$ and patch to an intermediate 
layer of material which respects energy conditions. This intermediatelayer 
may then be patched to the Schwarzschild-(anti) de Sitter vacuum and thus 
there will be little or no energy condition violation. Notice from figure 3 
that averaged energy conditions (energy conditions integrated over the stellar 
bulk) are certainly satisfied as most of the stellar bulk is "well behaved" in 
the context of energy conditions. Also interesting is that, at the throat, the 
matter is still exotic in the sense that the magnitude of the pressure is of 
the same order of magnitude as the energy density. A positive cosmological 
constant, if present, would minimize the amount of material needed at the 
throat although the cosmological constant would most likely be quite small 
when compared to 1/Yq and therefore would not necessarily eliminate this 
exotic behaviour. 

The weak energy condition may also most easily be studied by form- 
ing the quantity [|J, |§ 

e := ^+A. (48) 

\p\ 

This quantity should be positive for any matter respecting the weak energy 
condition. We plot this quantity for the example above in figure 3f. This 
figure shows that throughout the stellar bulk, with the exception of a small 
region near (but not at) the throat, the weak energy condition is satisfied. 
The violation, however, is very small. This agrees with the perturbative anal- 
ysis above, as well as with ||23|| , where it si shown that the energy violations 



can be confined to an arbitrarily small region around the throat provided 
that &i(r) is close to unity. As pointed out in much of the literature on 
energy conditions, this violation does not conflict with experimental results 
as the weak energy violating Casimir effect has been verified experimentally 



0- 



Another interesting example, which we will briefly consider, is that 
where the wormhole connects not just two universes but any number of uni- 
verses. An in-depth analysis of this situation will not be made and this 
example simply serves to illustrate that a minor extension to the above tech- 
niques is all that is required for such a system. 
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The profile curve of the "top half" of a single closed universe is shown 
in figure 4a. The bottom half of the universe may be generated by reflection 
via P(r) — ► —P(r). We use as an example of such a profile curve the following 
function: 



P(r) = — cos x 
k 



B) 



R 



(49) 



where k, R and B are constants (with B > R) related to the size of the 
universe and to the size of the throat via d = B+R and r = B—R << B+R. 
This profile, along with its reflection about the r-axis serves to generate the 
closed universe. Note that this universe may be joined at its throats to other 
(potentially infinite number of) closed universes or to an open universe at 
each end of the closed universe chain. In this latter case the closed system 
represents a set of "baby universes" to the open ones. It should also be 
possible to terminate the chain not with open universes but instead a closed 
universe with only one throat. Note that continuity of the metric must 
be respected at each coordinate chart junction in order to have a smooth 
patching. A representative system of three closed universes with similar 
radii is shown in figure 4b. 

Again we utilize the mixed method of solving the Einstein field equa- 
tions in order to minimize energy condition violation. Singularities and event 
horizons are still forbidden and, since here we do not patch to a well behaved 
vacuum solution, a slight modification of the above techniques will be re- 
quired to ensure these conditions both at the throat and at r = d :— (B + R). 

The Einstein equations give, for the energy density of the fluid, 

8vrp + A = —b{r) tl 

k 2 (R 2 + r 2 -B 2 ) + l 
r 2 [k 2 ((B - r) 2 - R 2 ) - if 

where the ± subscript has been dropped as it is implied that quantities in 
the top half of the universe should be equal to those in the bottom half. 
Although this condition is not a strict requirement of the model, it is used 
here to simplify the system. It should be noted that, as mentioned above, 
the metric must be continuous at both at the throat and at r = d in order 
to have a smooth patch. Also, the junction conditions require T\(d) in the 
lower half to equal T\ (d) in the upper half. 
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4b 



4a 




Figure 4: a) The profile curve for the upper section of a single universe, b) A 
chain of three closed universes connected by wormholes. 



Equation ( p0|) may be positive or negative. If one wishes to demand 
a positive total energy density everywhere, then the condition r < ^r^ must 
hold. However, the presence of a negative cosmological constant allows ( |50| ) 
to be negative and still have everywhere non-negative energy density for the 
fluid. 

We now analyze the manifold for event horizons and singularities. In 
this example, the Einstein equations give: 



7( r ),i 



T(r) + 



1 + 



k 2 [R 2 - (r - B) 2 ] 



(51) 



where the previous notation is still used (r(r) := 8ttT\ — A). Note that now 
there are two potentially troublesome spots, namely r = r$ and r = d. As 
before, to guarantee that the integral on ( [5lD is well behaved, it is required 
that r + 1/r 2 vanish at least as fast as the denominator of the second term 
in square brackets at these spots. Again, as is common with static wormhole 
solutions, a net tension is required at the throat to maintain the wormhole. 
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A tension is also required at r = d. Although many functions will possess 
this behaviour, for the specific example considered here we assume 

1 / r — ro inn \ 1 / d — r run 



r(r) = -^ cos Kd—7,— ) ~ ffi cos [d-7 — ) ' (52) 

where m is an odd integer. 

The Riemann curvature tensor in the orthonormal frame at r = r$ 
gives, from (|30a|) : 



R (trtr)(r ) = 4r 3 (ro _ 6(ro)) {r 2 (r )r 6 + r(r )rg (1 + 6(r„),i) 

+ r 2 fo(ro),i}- (53) 

By noting that b(r ) = B — R and 6(r ),i = 2k 2 R(R — B) + 1 it is easy to 
check that (B3) sets the quadratic in the numerator of the above expression 
to zero. Similarly, at r = d 

R(trt r) (d) = M3 (rf X _ {r 2 (d)r 6 + r(d)d 4 (1 + 6(d),i) 

+ d 2 6(rf),i}- (54) 

Here 6(d) = B + R and 6(d) ,i = 2k 2 R(R + B) + 1 and again (|2|) gives a 
finite result. Finally, it should be noted that, since p, r(r) and 7(r) : i are 
everywhere finite, the transverse pressures (defined via equation (|j)) are also 
well behaved. 



5 Conclusion 

In this paper a general model of a static wormhole system was developed. 
The construction was designed to encompass a large class of static wormhole 
solutions and limitations were imposed on it so that physical quantities are 
as reasonable as possible. The system was shown to obey both energy con- 
ditions at the throat as well as averaged energy conditions throughout the 
matter bulk. Also, the matter system smoothly joins a vacuum possessing an 
arbitrary cosmo logical constant. Many specific examples may be constructed 
from the work here and a couple were considered, including a system of mul- 
tiple connected universes. It would be interesting to derive a general class 
of time-dependent wormhole geometries, although the equations governing 
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such a system are formidable. In this case one would also need to specify, 
at the very least, some qualitative features of the future evolution of the 
system. Examples are cases where the wormhole effectively "closes up" or 
simply translates through spacetime. As well there is the issue of causality 
( [[25)1 , f26H , f27j, dm and others). There is also much interesting work to be 



done regarding wormholes of other symmetries (for example see [29], [[3C]). 
It should be straight-forward to apply the analysis here to generate a general 
class of cylindrical wormholes, for example. 

It is hoped that this work may provide a general wormhole model 
from which one may base future work in this fascinating field. The aim is 
also to shed light on methods for satisfying all field equations and identities 
in a way relevant to studies of stellar systems. 
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